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The complex nature of the nuclear star cluster in FCC 277* 
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ABSTRACT 

Recent observations have shown that compact nuclear star clusters (NSCs) are present in up 
to 80% of galaxies. However, detailed studies of their dynamical and chemical properties are 
confined mainly to spiral galaxy hosts, where they are more easily observed. In this paper we 
present our study of the NSC in FCC 277, a nucleated elliptical galaxy in the Fornax cluster. 
We use a combination of adaptive optics assisted near-infrared integral field spectroscopy, 
Hubble Space Telescope imaging, and literature long slit data. We show that while the NSC 
does not appear to rotate within our detection limit of ^6 km s~\ rotation is detected at larger 
radii, where the isophotes appear to be disky, suggesting the presence of a nuclear disk. We 
also observe a distinct central velocity dispersion drop that is indicative of a dynamically cold 
rotating sub-system. Following the results of orbit-based dynamical modelling, co-rotating as 
well as counter-rotating stellar orbits are simultaneously needed to reproduce the observed 
kinematics. We find evidence for varying stellar populations, with the NSC and nuclear disk 
hosting younger and more metal rich stars than the main body of the galaxy. We argue that gas 
dissipation and some level of merging have likely played an important role in the formation of 
the nucleus of this intermediate-mass galaxy. This is in contrast to NSCs in low-mass early- 
type galaxies, which may have been formed primarily through the infall of star clusters. 

Key words: Galaxies: elliptical and lenticular, cD; Galaxies: formation; Galaxies: nuclei; 
Galaxies: kinematics and dynamics 



1 INTRODUCTION 

It is now believed that up to 80% of all galaxi es host nuc- 
lear s tar clusters (NSC s ) in their centr es (e.g. ICarollo et al.l 
1 19981 : iBoker et all 120021 : ICote et alj l2006h . Typically, low- and 
intermediate-luminosity galaxies show a central light excess inside 
a characteristic radius rb~ 0. 02_R e q above t he inner extrapola- 
tion of the global light profile dCote et alj|2007h . Th ese NSCs usu- 
ally reside in the photom etric centre of the galaxy dBinggeli et al.l 
l2000l : lB6keretal.ll2002h and their loc ation overlaps with the kin- 
ematic centre jNeumaver et al.l 1201 it) . NSCs are usually brighter 
than typical globular clusters, compact (r ~ 5 pc), massive 
(M ~ 10 7 Mq), may be flattened, and often contain multiple stel- 
lar populations and complex structures (e .g. Walcher et alj 200a 



lar populations and complex structures (e.g. Walcher et al. 2UUD, 
l2006t ICote et alj|200r3 : iRossa et alj|2006l : ISeth et alj|200q booa 
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iBarth et alj|2009l ; iTurner et alj 120121: iPiqueras Lopez et al.ll2012h . 
The ir masses seem to co r relate with the mass o f the host galax- 
ies dFerrarese et alJlioOcjt IWehner &~H arris 2006), extending the 
super-massive black holes scaling relations to the low-mass end 
of galaxies. In some cases , NSCs appear to co-ex ist with central 
black holes (e.g. revi ew bvlGraham & Spitler 2009, and references 
therein) and recently iNeumaver & Walchen (120 1 2I> suggested that 
NSCs may be the precursors of massive black holes in galaxy nuc- 
lei. However, only handful of detailed studies on the properties of 
NSCs exists and these are mainly focused on such objects in late- 
type galaxies. Characterising NSCs in early-type galaxies is a non 
trivial task, both because of the high surface brightness of the un- 
derlying galaxy, and because of the NSCs compact sizes. 

With the availability of adaptive optics fed integral field unit 
(IFU) instruments, this task is now becoming feasible. For ex- 
ample, ISeth et af have shown that the nucleus of NGC 404, 
a nearby SO galaxy, hosts several morphologically and dynamic- 
ally distinct components. The NSC in this galaxy shows a modest 
rotation aligned with the galaxy, a gas disk that rotates perpendicu- 
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Figure 1. HST/ACS z-band image of FCC 277. The star that we used as a 
NGS for the AO correction is masked out here. The SINFONI FoV (3" X 
3") is indicated with a white box. The overlaid contours are from the MGE 
light model, discussed in Sect. [5] 



larly to the stars, and probably an intermediate-mass black hole (~ 
10 s Mq). Such complicated structure inevitably poses the question 
of how NSCs have formed. Currently, there are two main scenarios 
proposed. The first involves the dissipationless infal l of star clusters 
to th e galaxy centre due to dynamical friction dTremaine et al.l 
Il975t) . The second suggests NSCs t o be the result of dissipat ional 
sinking of gas to the galactic centre dMihos & Hernquisal 19941) . 

Numerical simulations of globular clusters infall have had 
certain success in reproducing the observed surface bright- 
ness profiles of nucleated galaxies , although with larger nuc- 
lei sizes comparing to what is obs erved (e.g. lOh & Lin] |2000| ; 
Capuzzo-Dolcett a & Miocchi 2008a b). However, more recently 
Hartmann et al. 42011b showed that star cluster accretion onto a 
pre-existing nuclear disk did not produce the observed line-of- 
sight kinematics of NSCs. They suggested that purely stellar dy- 
namical mergers cannot be solely responsible for the formation 
of NSCs and that gas dissipation must also play a significant 
role in assembling the cluster's mass. What is the exact origin 
of this gas and h ow it gets transpor ted to the galaxy nucleus is 
still under debate. iBekki et al] J2OO6) have shown that the dissip- 
ative merging of stellar and gaseous clumps formed from nuclear 
gaseous spiral arms in a gas disk eventually produce nuclei that 
rotate, are flattened and have a ran ge of ages and metallicities. 
Pfl amm-Altenburg & Kroupal d2009h concluded that compact star 
clusters with masses ^ 1O 6 M0 act as cloud condensation nuclei 
and are able to accrete gas recurrently from a warm interstellar me- 
dium. This may cause further star formation events and account for 
multiple stellar populations in the most massiv e globular and nuc- 
lear star clusters. Recently. iTurner et al.l (120121) concluded that the 
dominant mechanism for nucleus growth in low mass early-type 
galaxies is probably infall of star clusters through dynamical fric- 
tion, while at higher masses, gas accretion resulting from mergers 
and torques becomes dominant. 

In this paper we present a detailed study of the nucleus in the 
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Figure 2. An HST/ACS z-band profile (left panel) and (g~z) colour (right) 
of FCC 277 iTurner et al.l201 2; Ferrarese et al 12013b . The coloured curves 
show the Sersic fits to the two components: nucleus and galaxy. The arrows 
at (f.'25 indicate the break radius, r^, at which the nucleus component starts 
to dominate the galaxy light. 

intermediate-mass early-type galaxy FCC 277 (NGC 1428) and dis- 
cuss our observations in the light of the current assumptions of nuc- 
leus formation. This galaxy is a member of the Fornax cluster and 
is part of the ACS Fornax Cluster Survey (Jorda n etal.ll2007h . Its 
basic properties, as well as the main parameter of the NSC, are 
listed in Tabled In Fig. □ we show part of the HST/ACS 2-band 
image, together with the field-of-view of VLT/SINFONI that we 
used to complete our study. In Fig. [2] (left panel) we plotted the 
z-band surface brightness profile, together with the two Sersic fits 
(dashed lines) that describe the galaxy light. The outer galaxy light 
is represented with a Sersic fit with n = 1.8. The point where the 
nucleus starts to dominate over the inner extrapolation of the Sersic 
fit is called break radius and for FCC 277 has the value r^= 0'.'25 
(indicated with an arrow in Fig. [2}. The nucleus is fitted with an- 
other Sersic profil e with n = 1 . 7. Fo r a detailed descr i ption of the 
fitting process see lTurner et al. l d2012l) ; lFerrarese et al.l (120131) . 

The paper is organised as follows: in Sect. |2]we describe our 
observations and data reduction. In Sect. [3] we compare the light 
profiles obtained from HST/ACS and VLT/SINFONI coupled with 
adaptive optics. Sect.|4]is devoted to the kinematics analysis of our 
IFU data, while in Sect. [5] we present results from dynamical mod- 
elling of the galaxy and the NSC. In Sect.[6]we explore the stellar 
populations of the nucleus of this galaxy. In Sect. [7] we discuss our 
findings in the light of current galaxy and NSC formation models. 
We conclude in Sect. [8] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

We obtained integral field spectroscopy of the Fornax E5 galaxy 
FCC 2 77 (NGC 1428) u sing VLT/SINFONI teisenhauer etal] 
120031 : iBonnet et alj[2004h in Natural Guide Star adaptive optics 
mode on October 6, 7, and 10, 2007 (programme ID 380.B-0530, 
PI L. Infante). We used the if-band grating (1.95 - 2.45 /im) that 
gives a spectral resolution R-~3500 (6.2 A FWHM as measured on 
sky lines). Our observations cover the central 3"x 3", with a spa- 
tial sampling of (/.'05 x O'.'lO. As a natural guide star we used a 
R = 14 m star located at 3'.'5 to the North of the galaxy centre (see 
Fig.QJ. Due to its proximity to the galaxy centre, this star does not 
appear in the guide star catalogue or the USNO catalogue as a sep- 
arate entry. Its celestial coordinates are a(J2000) = 03:42:22.9 
and r5( J2000) = -35:09:11.4. Our observations were carried out in 
service mode. 
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Table 1. Basic properties of FCC 277 and its Nuclear Star Cluster. 



FCC 277 




Reference 


Morphological Type 


E5 


Fereuson (1989) 




13. 8 m 


// 


Effective radius 


10.2" 


// 


(g — z) colour 


1.31±0.01 


Blakeslee et al. (2009) 


Distance 


20.7±0.7 Mpc 


// 


Major axis position angle 


115° 


Graham et al. (1998) 


Velocity dispersion 


81.7 km s" 1 


Weener et al. (2003) 


M vir 


~ 8 x 1O 9 M 


l) 


NSC 


Effective radius 


0.09" ~ 9 pc 


Turner et al. (2012) 


9 (mag) 


20.08 m ±0.16 


// 


(»-*) 


1.33±0.18 


// 



) Using M vir = 5.0 R cf! a 2 /G jCappellari et al.l200rj) 
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Figure 3. Comparison between isophotal parameters obtained from space 
and adaptive optics ground-based imaging. The red symbols show the light 
profile of FCC 277 as derived from our SINFONI A'-band reconstructed 
image. The black symbols denote HST/ACS z-band imaging. The two 
vertical dashed lines in panel (a) show the area where the two profiles were 
normalised (0"5 < r < C/7). The vertical arrows in each panel indicate the 
break radius, as in Fig. [2] In panel (b) we compare the 4 th cosine term of 
the isophotal fits that is indicative for deviations from pure elliptical shape. 
In panel (c) the measured position angles are displayed, and in panel (d) the 
ellipticity of the isophotes is shown. 

For the observations we used the standard near-IR nodding 
technique. Each observing block consisted of a sequence of object 
and sky frames (OOSOOSOOS), each individual integration was 
300 s, the sky fields were offset by 50" to the North. Science frames 
were dithered by 0'.'05 and 07 15 in order to reject bad pixels. There 
were six observing blocks. The total on-source integration time was 
3 hours. Additionally, after each observing block and at a similar 
airmass, we observed a B dwarf to act as a telluric star. 

2.2 Data reduction 

We used the ESO SINFONI pipeline v2.0.5 to perform the basic 
data reduction on each observing block, consisting of six object 
and three sky exposures. In brief, the pipeline extracts the raw data, 



applies distortion, bad pixels and flat-field corrections, wavelength 
calibration, and stores the combined sky-subtracted spectra from 
one observing block in a 3-dimensional data cube. For each resu lt- 
ing data cube, we then ran the lac 3d code toavies et al.ll201oT) to 
detect and correct residual bad pixels identified using a 3D Lapla- 
cian edge detection method. 

We reduced the telluric stars in the same way as the science 
frames. Then for each telluric star we extracted a one-dimensional 
spectrum, removed the hydrogen Brackett7 absorption line at 
2.166 /im after fitting it with a Lorentzian profile, and divided the 
star spectrum by a black body spectrum with the same temperat- 
ure as the star. The last step in preparing the telluric spectrum was 
to apply small shifts (<0.05 pixels) and scalings to minimise the 
residuals of the telluric features. To do this, we extracted a central 
one-dimensional spectrum from each science data cube and cross- 
correlated and fitted it with the corresponding telluric spectrum. 
Then we divided each individual spaxel in the six galaxy data cubes 
by the corresponding best fitting telluric spectrum. In this way we 
also obtained a relative flux calibration. 

Finally, we combined the six galaxy data cubes, using a 
3cr-clipping pixel reject algorithm. We also reconstructed a two- 
dimensional image of the galaxy, after integrating the spectral di- 
mension of the final data cube in the range 2.1 - 2.4 /im, where 
the contamination from sky lines residuals is minimal. To be able 
to robustly measure the velocity and velocity dispersion from the 
spectra of the galaxy a minimum signal to noise of 20 per pixel 
is required. Because the SINFONI pipeline does not provide error 
propagation during data reduction, we estimated the noise in each 
spectrum of the data cube as the r.m.s. of the residuals after sub- 
tracting a smoothed model of the spectrum. Then we used this noise 
estimate to bin the final galaxy data cube to achieve an approxim- 
ately constant S/N~25 us ing the Voronoi 2D binning method of 
ICappellari & Copinl 12003). This S/N allowed us to conserve a good 
spatial resolution (our smallest bins in the centre are ~0'.' 1 across), 
while we are still able to reliably extract the stellar kinematics. 



3 COMPARISON BETWEEN VLT/SINFONI AND 
HST/ACS LIGHT PROFILES 

In Fig. [3] we compare the light profiles of FCC 277 derived from 
HST/ACS imaging and VLT/SINFONI ob servations. T his galaxy 
is part of the ACS Fornax Cluster Survey (Jordan et al. 2007|) and 
high spatial resolution imaging in the g-band and z-band from the 
ACS is available. We used the IRAF task ellipse to fit elliptical iso- 
photes to the ACS z-band image (black symbols in Fig. [3) and to 
the SINFONI reconstructed image (before binning; red symbols). 
In panel (a) the two luminosity profiles are compared, after being 
normalised in the region 0'.'5 - 077 (dashed vertical lines). Isophotal 
parameters are plotted against the semi-major axis length. The ver- 
tical arrows denote the break radius, 0'.'25, at which point the 
nuclear compone nt starts to dominate the surface brightness profile 
dCote et alj200?t see also Fig.|2]l. 

Although observed with adaptive optics, the SINFONI light 
profile is less steep than the ACS profile within the inner ~0'.'3. 
Assuming this is due to the lower spatial resolution of the SINFONI 
data, we estimated their PSF to be 0'.' 16 5 (FWHM) by convolving 
the ACS image (using a Tiny Tim PSF. iKrisnl 1995T) with a given 
Gaussian PSF until it matches the light distribution in the SINFONI 
image (observations of stellar PSFs were not obtained during the 
SINFONI run). 

In panel (b) of Fig. [3] we plotted the cosine 4 th order Fourier 
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Figure 4. The spectrum from the central bin (S/N~25) of FCC 277 with an 
over-plotted best fitting template spectrum, derived by the pPXF code (red 
line; fit residuals are shown in grey). The location of the strongest near-IR 
absorption features in the X-band is indicated. The dashed vertical lines 
show the location of the strongest sky emission lines. 
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Figure 5. Velocity (top panel) and velocity dispersion (bottom panel) maps 
of FCC 277 and corresponding errors (right panels). Over-plotted are con- 
tours with constant surface brightness, as derived from our reconstructed 
SINFONI if -band image. 



coefficient of the isophotes, divided by the semi-major axis length. 
Positive values of this parameter are indicative of disky isophotes, 
as are observed in both the z- and 7\"-band images in the inner 1". 
There are two peaks in the a4/a profile, one at ~(X'2, coinciding 
with the break radius, and a stronger second peak at ~0'.'6, coin- 
ciding with the peak in the velocity field (see Sect.[4]l. In panel (c) 
no significant variations of the position angle outside of the break 
radius are observed and the mean PA is consistent with the one de- 
rived at larger radii (see Table [T]. In panel (d) the ellipticity reaches 
a maximum at ~0'.'35 for both the ACS and SINFONI profiles, al- 
though with different amplitudes. These differences are expected 
due to the differences in t he PSF of the two images; a larger PSF 
leads to rounder isophotes dPeletieretal.ll990l) . The comparison of 
the two profiles led us to the conclusion that the SINFONI ground 
based adaptive optics assisted observations are similar in quality to 
the HST/ACS images. 

The observed features of the isophotal parameters point to a 
picture where, within the break radius, the nuclear star cluster may 
be flattened or, alternatively, may be the superposition of a round 
NSC and a larger scale disk. Such nuclear disk beyond the break 
radius is evident i n the diskiness parameter aA/a at ~0'.'6 (see also 
iTurner et al.l2012h . 



4 STELLAR KINEMATICS 



We used the pPXF code JCappellari & Emsellemll2004l) to derive 
the first and second order of the line-of-sight velocity distribution, 
working with a library of seven template spectra of K and M gi- 
ant stars. These templates were observed with the same instrument 
and the same setup as our science target. To find the best fitting 
composite template spectrum we used the region between 2. 1 and 
2.36 /im, where several strong absorption features allow accurate 
measurements (see Fig . [4}, and masked the strong near-IR sky lines 
dRousselot et alj2000f) . In Fig.|4]we show the spectrum of a central 
bin in the galaxy with the over-plotted best fitting composite tem- 
plate as derived by the pPXF code (in red), as well as the residuals 
(in grey). 

Our stellar mean velocity and velocity dispersion maps are 



40 
20 

-20 
-40 
120 



100 




5s _*_•[£ * 



• This work 
A Koleva et al. 201 1 
.□ Spolaoretal. 2010 
O Graham etal. 1998 

— I 1 

. NW 






-10.0 



-1.0 



-0.1 0.1 
rad [arcsec] 



1.0 



10.0 



Figure 6. Behaviour of the ste llar kinematics at different spatial scales. 
Long-slit kin e matics data from Graham et al. ( 1998, diam ond symbols), 
Spolao ret al.1 d2010l blue squares) and iKoleva et alj 1201 ll red triangles) 
are shown. With solid symbols we plotted the innermost kinematics profile 
of FCC 277, extracted from the SINFONI maps using kinemetry. 



shown in Fig. \5\ Using kinemetry, described by iKrainovic et al.l 
(2006), we extracted the velocity and velocity dispersion profiles, 
shown with filled symbols in Fig. [6] We observe rotation around 
the minor axis of the galaxy up to ±25 km s~ at r ~0 / .'6, which 
is outside the break radius of the luminosity profile (rb= 0'.'25, 
see Fig. [2] left panel) where the NSC is supposed to reside. 
At the centre of the galaxy the velocity dispersion approaches 
~55 km s _1 and then it increases in the outer parts of the field of 
view to reach ~90 km s _1 at r ~ 1". 

In Fig. [6] we also compare our own SINFONI data with a 
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compilation of literature measurements, der ived using long slits , 
aligned along the major axis of the galaxy. iGraham et al.l d 1 998L 
diamond symbols) used a spectrograph with a slit width 2" on the 
Australian National Universitys 2.3 m telescope at Siding Spring 
Observatory. Their data do not cover the inner 4 "of the galaxy, due 
to the presence of a relatively bright star clos e to the nucleus (that 
we used as a natural guide star for the AO). ISpolaor et all d20ld 
blue squares) used the GEMINI/GMOS instrument with a slit width 
of 1". The seeing durin g these observations was in the range 0'.'7 - 
l". lKoleva et alj (201 ll . red triangles) reanalysed the same observa- 
tions. 

Based on Fig. [6] we conclude that the rotating substructure in 
the nucleus of FCC 277 co-rotates with the main body of the galaxy 
and that the velocity dispersion in the outer regions of the SIN- 
FONI field-of-view reaches similar values as the long slit studies. 
We note that the spatial resolution of our adaptive optics supported 
data is much higher than the resolution achieved by the other three 
studies, thus we cannot directly compare the radial profiles at small 
galactocentric radii. The significantly worse spatial resolution of 
the long-slit observations means that the inner rotation and velocity 
dispersion dip are washed out. 

The observed rotation, taken together with the drop in the ve- 
locity dispersion, indicates the presence of a co-rotating cold sub- 
structure in the inner 0'.'6 of the galaxy. We fitted the kinematic 
position angle of this substructure usin g the method described 
in Appendix C of iKrainovic et alj d2006l) . The measured value is 
11 8° ±10°, which is consistent with the photometric position angle 
of the main body of the galaxy, derived at larger radii (see Table[T](. 

For early-type galaxies the apparent stellar angular mo- 
mentum A_r c and the galaxy flattening are now a commonly 
used tool to classify galaxies into fast- and i/ow-rotators 
( lEmsellem et"alll2007l : Cappellari eud"1l2007h . The method needs 
IFU data to measure \r c inside one-effective radius. For FCC 277 
only long slit data are available out to the effective radius. Thus 
we used our best-fit Schwarzschild model from Sect.[5]to simulate 
the velocity and velocity dispersion as they would be observed by 
an IFU. We measured Xr, = 0.3. Using the most recent classi- 
fication from lEmsellem et all ( 1201 it) , we found that this galaxy is 
a/asf-rotator and lies slightly above the dividing line between the 
two classes, with its ellipticity e = 0.3. 



5 DYNAMICAL MODELLING 

To measure the mass distribution and the orbit configu ration of the 
inner part of FCC 277 we used Schwarzschildl d 19791) modelling. 
This method ( Ivan den Bosch et alj 20081) works by constructing a 
trial mass model of the galaxy, including a black hole, stars and 
dark halo. Then, the gravitational potential is inferred from the 
mass model and representative orbits are integrated numerically, 
while keeping track of the paths and orbital velocities of each or- 
bit. We can then create a mass model of the galaxy by assigning an 
amount of mass to each orbit so that the overall stellar mass dis- 
tribution is reproduced, while simultaneously fitting the observed 
stellar kinematics. The effect of the PSF on the observed stellar 
kinematics is an integral part of the dynamical model. These mod- 
els have the advantage that they do not require any assumptions 
regarding the orbital anisotropy of the galaxy. 

The models were constructed as follows. First, we paramet- 
rized the galaxy stellar surface brightnes s using the multi-G aussian 
expansion (MGE) method, described bv lCappellaril d2002l) , on the 
ACS z-band image. In Fig. Q] we show this image with overlaid 




Figure 8. Confidence interval of the dynamical models of FCC 277 for the 
mass-to-light ratio in z-band. The black dots indicate the location of the 
models and the contours indicate 1, 2, and 3cr intervals, where the 3-sigma 
level is indicated by a thick line. 



contours of the MGE light model. There were 13 Gaussians with 
varying flattening fitted in total, the first two of them describing 
the NSC. The galaxy shows strong rotation around the minor axis 
and we therefore assumed the galaxy is obla te axisymmetric , whic h 
is the most common configuration (e.g. [Padilla & St rauss! Eo08l) . 
The galaxy is also strongly flattened, with a minimum flattening 
0.6 at 15", and can thus not be seen more face-on than i = 65°. 
Then we used our symmetrised (using the meth od described in 
Appendix A of Ivan den Bosch & de Zeeuwll2010l) SINFONI IFU 
kinemat ics from Sect.[4]and the Sc hwarzschild orbit superposition 
method ( Ivan den Bosch et ap2008l) to construct a realistic dynam- 
ical model for the galaxy and th e NSC. We also included the long- 
slit data of Graham et al.1 dl998l) . to be able to constrain the mass- 
to-light ratio of the main bod y of the galaxy. We did n ot include the 
data of ISpolaor et al.1 J2010b and lKoleva et alj J201lh . which result 
from two different reductions of the same data set, because they do 
not match with our kinematics measurements for the inner parts of 
the galaxy. We probed the following parameters: the central black 
hole mass and separate mass-to-light ratio (M/L) for the galaxy and 
the NSC, using 5000 dynamical models. Changes in the inclination 
between 65°and 90° led to insignificant changes in the M/L, thus 
we marginalised over it. 

In Fig.|7]we show the input symmetrised SINFONI kinematics 
velocity and velocity dispersion maps (left panels) together with the 
resulting kinematics obtained by the best fitting dynamical model 
(middle panels). The reduced \ 2 of the best models is ~0.21 over 
the 81 SINFONI bins. The low value of the reduced \ 2 ls due to 
the very conservative estimate of our kinematics errors. The best-fit 
M/L z of the galaxy and NSC is 3.2 ± 0.4 and 3.0 ± 1.0 respect- 
ively, as shown in Fig. [8] Confidence intervals are determined using 
Ax 2 statistics, assuming two degrees-of-freedom. Thus, the mass 
of the NSC is 1.4 ± 0.4 x 10 7 M Q . The black hole mass is uncon- 
strained, as the uncertainties on the central kinematics are too large. 
The difference in velocity dispersion between a 10 5 and a 1O 7 M0 
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Figure 7. Comparison between our symmetrised SINFONI kinematics (left panels) and the ones obtained by the best fitting dynamical model (middle panels). 
The right panels show the resultant kinematics maps from our best fitting model when we do not include stars on counter-rotating orbits. 



black hole is 5 km s _1 and the uncertainties on a are ~ 15 km s -1 . 
Black hole masses above 1O 7 M0 do yield significantly worse fits 
and hence we place an upper limit of 10 7 A/q. To robustly determ- 
ine the black hole mass higher S/N spectra of the nucleus need to 
be obtained to reduce the uncertainties on the central kinematics 
and the PSF of the IFU data needs to be known precisely. The in- 
clusion of a dark matter halo does not alter the M/L of the NSC 
and only very weakly the M/L of the galaxy. The M/L m \ is ex- 
pected to contain only a small contribution from the dark matter 
(Cappellari et al. 2006), hence our final adopted model does not in- 
clude a dark matter halo. Kinematics reaching much further out are 
need to properly constrain the dark matter halo. 

Apart from the mass distribution of the galaxy, the models 
yield the orbital distribution as a function of radius. In Fig. [9] we 
show the orbital mass weights as a function of the average radius 
and spin \ z — J z x (f/ff), where J z is the average angular mo- 
mentum along the short z-axis and a the average second moment 
of the orbits. 

We detect the presence of three distinct components: both a 
co- and counter-rotating component as well as a non-rotating bulge 
component. The relative contribution of each of these components 
is shown in the bottom panel of Fig. [9] Both rotating components 
extend well inside the break radius and have similar contributions 
in the NSC region. The sigma drop seen in the stellar dispersion 
map coincides with a decrease of the non-rotating orbits. The ques- 
tion arises if this is the only possible orbital configuration for this 
system. The orbits available to the Schwarzschild models are a fully 
representative set and the linear solver used to construct the models 
is guaranteed to find the g lobal minimum ( Ivan de Ven et al ] |2008l ; 
Ivan den Bosch et al]|2008l) . This guarantees that the model finds 



the best-fitting orbital configuration. There could be other solutions 
that are also a good representation of the observations. As a con- 
sistency check, we attempted to fit models without counter-rotating 
orbits, which would exclude an opposite angular momentum accre- 
tion event as a formation scenario for the nucleus. However this led 
to a significantly worse match of the stellar kinematics (see Fig. [7] 
right panels), which indicates that counter-rotating orbits are thus 
required. 

We note that we were unable to use the Jeans' modelling ap- 
proach to fit the stellar kinematics using the method described by 
ICappellaril 120081) . Although we could receive reasonable fits to the 
second velocity moment (Vrms), the fits to the velocity field were 
unsatisfactory. This is because, as it is at the moment, the JAM 
package does not allow the rotation parameter k to accept positive 
and negative values simultaneously for a given MGE Gaussian. 

To be able to quantitatively discuss the different orbital frac- 
tions, the dark matter halo of the galaxy and its global M/L, one 
would need improved long-slit or other large scale kinematics. 



6 STELLAR POPULATION PARAMETERS 

In addition to the structure and dynamics, we can investigate what 
the parameters of the stellar populations of the NSC and/or nuclear 
disk in the heart of FCC 277 are and whether they differ from the 
main body of the galaxy. 

Usually, nuclei in low luminos ity Fornax and Virgo galaxies 
are b luer compared to their hosts (ICote et al .1120061 iTurner et al.l 
l2012h . In the right panel of Fig. [2] we show the (q — z) colour 
profile, as derived from HST/ACS imaging l lTurner et al.l 120121 
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Figure 9. Top panel: Distribution of mass along the orbits of our best fitting 
dynamical model, as a function of angular momentum and radius. Colour 
coding reflects a factor of 3.5 span in mass density (darker colour corres- 
ponds to higher mass). Bottom panel: Fraction of non-/rotating orbits as 
a function of radius. A non rotating bulge (—0.1 < A z < 0.1) is de- 
noted with a black line. The red line corresponds to the sum of all co- 
rotating orbits with X z > 0.1, the blue line - all counter-rotating orbits 
with X z < -0.1 



iFerrarese et al.l 12013b . The integrated colour of the nucleus is 
(g - z) = 1.33 ± 0.18 dTurneretal.ll2012l) and does not dif- 
fer from the main body of the galaxy. If nuclei follow the same 
colour -metallicity relat ion as globular clusters in early -type galax- 
ies do JPeng et alj2 006), than the red colour would be indicative for 
higher metallicity of the NSC. However, age effects cannot be ex- 
cluded, due to the well known age-metallicity degeneracy of broad 
ban d colours. 

iKoleva et alj 1201 ll) measured from optical spectroscopy the 
age and metallicity of the core of FCC 277 (within a 0>.'5 radius 
aperture) to be 5.4 Gyr and [Fe/H]=— 0.07, respectively. At the ef- 
fective radius these values are 7.7 Gyr and [Fe/H]=— 0.50. Their 
data lack the spatial resolution to differentiate the NSC and the disk, 
however there is a pronounced negative age and positive metallicity 
gradient towards the nucleus. 

We used our near-IR IFU spectra to measure the line strengths 
of Nal (~2,2 ^m) and 12 CO(2-0) (~2.3 /itn) absorption features 
(see Fig.|4j- From previous stellar population studies in the near-IR 
wavelength range we kn ow that the Nal index increases with metal- 
licity and younger age dSilva et al. 2008; Mar mol-Oueralto et al.l 
120091 : iLvubenova et aljfeoid) . and the Deo index is e xpected to 
increa se with higher metallicity for ages above 3 Gyr dMarastonl 
120051) . 

We measured the two indices us ing the definition of 
lFrogeletal.l ( l200ll) for the Nal index and Marmol-Oueral to et al.l 
d2008l) for the Deo index. Before measuring, we first broadened 
our spectra to 6.9 A (FWHM, ~ 94 km s~ x ) to match the spectral 
resolution of othe r stellar population studies of elliptical galaxie s 
in the near-IR (e.g. lSilva et alj200ctlMarmol-Oueralt6 et alj2009h . 
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Figure 10. Nal and Deo index maps with HST like spatial resolution (left 
panels) and radial profiles {right panels). The red solid lines are the least 
squares linear fits to the data. 



Finally, we corrected the Nal index to ze ro velocity dispers ions us- 
ing the velocity dispersion corrections of lSilva et alj J2OO8). 

In Fig.[l0]we show our index maps as well as their radial pro- 
files. We observe radial gradients for both indices. Increase in Nal 
towards the centre is consistent with increasing metallicity and/or 
younger age in the nuclear disk and NSC. On the Deo map we 
see that the strongest index values towards the centre seem to form 
an elongated shape, aligned with the major axis of the galaxy and 
the rotating structure, visible on the velocity map (Fig. [5]) This in- 
crease is again consistent with the redder colour and indicative of 
higher metallicity. 

In Fig.QT]we compare FCC 277 with other early type galaxies 
in the Fornax cluster in terms of their Nal and D eo indices versus 
their central velocity dispersion. The spectra of ISilva et al.1 (2008) 
cover 1/8 of the effective radius of each galaxy and are marked with 
circles. Open circles represent galaxies with old stellar populations, 
solid symbols stand for the galaxies that have optical signatures 
of recent (<3 Gyr) star formation. The dashed line illustrates the 
least-squares linear fit to the old galaxies only. Our SINFONI ob- 
servations of FCC 277 cover approximately the same radial extent 
as the other galaxies. We plotted the values measured on the integ- 
rated spectrum with solid blue squares. We extracted the NSC area 
(r ^ 0'.'25) and marked our measurements with orange asterisks. 
With open red diamonds we indicated the values for the so called 
"nuclear disk", integrated over the range 0'.'25 < r 0'.'7. 

The Nal index of the NSC is much higher compared to the 
extrapolation (shown with a dotted line) of the a-Nal relation for 
old Fornax galaxies and is closer to the systems with younger 
(<3 Gyr) stellar populations. There is not a big difference of the 
index strength between the nuclear disk and the galaxy as a whole, 
so one would infer purely old age for these two components if look- 
ing only at this plot. The Deo index seems to saturate for galaxies 
with velocity dispersion higher than 100 km s -1 . Until reliable stel- 
lar population models for the near-IR become available, we cannot 
provide a quantitative estimate for the changes in stellar population 
parameters in the nuclei of early-type galaxies. At the moment we 
can only speculate that, taken together with the red colour of the 
NSC (as red as the host galaxy, which is unusual for NSCs in For- 
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Figure 11. Nal and Deo indices versus the central velocity dispersion of 
Fomax early-type galaxies. Open symbols denote old Fornax galaxies, solid 
symbols - ga laxies with optical signatures of recent (<3 Gyr) star formation 
I Silva et al. 2008). The dashed line represents the least-squares linear fit to 
the old Fornax galaxies only, the dotted line denotes its extrapolation to 
lower c. The NSC in FCC 277 is marked with an orange asterisk, the so 
called nuclear disk - with a red diamond, and the whole galaxy, as observed 
with SINFONI, is marked with a blue square. 

nax as discussed above), these point to a mixture between younger 
age and higher metallicity compared to the main body of the galaxy. 



7 FORMATION MECHANISM OF THE NSC 

So far we have collected evidence that: i) the nuclear star cluster 
inside the central 0'.'25 of FCC 277 does not rotate within our er- 
ror of ~ 6 km s" 1 ; if) the rotation visible at ~0'.'6 overlaps with 
a maximum in the diskiness of the isophotes, suggesting the exist- 
ence of a nuclear disk around the NSC; Hi) the existence of a central 
velocity dispersion drop is indicative of significant rotation in the 
same area; iv) dynamical modelling reveals a significant rotation 
in the inner 1" in both directions, i.e. co- and counter-rotation, ex- 
plaining the low level of observed rotation, the low \r c , and the 
sigma drop; v) there is no significant difference in the derived dy- 
namical M/L between the NSC and the galaxy. However, within 
the errors, a change in the stellar population parameters may not 
lead to an obvious change in the M/L for the different compon- 
ents. The derived dynamical M/L z within the errors is consist- 
ent with the predictions of stel lar population models for a Salpeter 
IMF and 5-10 Gyr age (e.g. lBruzual & Charlotll2003t iMarastonl 
2005). If one compares to Chabrier or Krpupa IMF, then our dy- 
namical M/L z is about a factor of two larger ; vi) there is evidence 
for differences in the stellar population parameters of the NSC, the 
nuclear disk, and the galaxy. This points to a scenario where the 
nucleus is younger and more metal rich. 

All these evidences for complex kinematics and stellar popu- 
lations point to a scenario where the NSC and disk formed through 
multiple episodes of gas accretion and subsequent episodes of star 



formation. Counter-roation points towards merge rs with orb i tal an- 
gular momentum opposite to the host galaxy. ISeth et aL I fcoioh 
reached the same conclusion for the nucleus in NGC 404, a nearby 
SO galaxy. Thus, FCC 277 is the second early-t ype galaxy that has a 
NSC exhibiting complex star formation history. iTurner et al 
noted that the lowest-mass galaxies with nuclei in the ACS Fornax 
and Virgo Cluster Surveys seemed to be structurally simple, having 
likely formed through star cluster infall. The more massive galaxies 
(such as FCC 277) seemed to be more structurally complex in their 
inner regions. They interpreted that as evidence for an increased 
importance of gas infall at higher masses. Our results about the nuc- 
leus of FCC 277 are consistent with this picture. Certainly, a larger 
sample is needed to study in more detail whether gas dissipation 
is a common mechanism for NSCs growth in ear l y-type galaxies, 
as it i s in late-type ones (e.g. IWalcher et al.ll2006l : iHartmann et al.l 
1201 lb . The NSC accounts for ~ 0.2% of the total mass of FCC 277, 
thus it is not an outlier of the typi cal scaling relations in early-type 
galaxies (e.g. lFerrarese et al1l2006l) . 



8 CONCLUDING REMARKS 

In this paper we present a pilot study of the detailed properties of 
nuclear star clusters (NSCs) in early-type galaxies. Although the 
nucleation frequency in galaxies is estimated to be ~ 80%, de- 
tailed data about the chemical and dynamical properties of NSCs 
exist mainly for spiral hosts. In early-type galaxies this task ob- 
servationally is not trivial due to the intrinsic brightness of the un- 
derlying galaxy light, as well as small angular extent of the NSCs, 
which at the distance of Fornax, span ~0 / .' 1 in diameter. We showed 
that using current technology, it is indeed possible and valuable 
information about the formation mechanisms of NSCs can be ob- 
tained. As a pilot test case we chose the galaxy FCC 277, a nucle- 
ated early-type galaxy that belongs to the Fornax cluster. This is the 
only member of this galaxy cluster that has a conveniently located 
bright star that one can use as a natural guide star for the adaptive 
optics (AO) system. 

Using SINFONI AO assisted observations, we observed the 
central 3" x 3" of the galaxy. Thus we obtained maps of the stel- 
lar kinematics with HST-like spatial resolution of 0'.' 165 (FWHM). 
Our velocity map (Fig. [5} shows clear rotation with a maximum at 
~f/.'6 from the galaxy centre, which overlaps with a maximum in 
the diskyness of the fitted isophots on the galaxy image (Fig. [3). 
The NSC itself has an effective radius of 0'.'08 and does not rotate 
within our detection limit of ~ 6 km s _1 . However, we observe a 
pronounced drop in the velocity dispersion in the central 1" that 
suggests the existence of a dynamically cold rotating sub-structure. 
Our dynamical modelling reveals that the nucleus of this galaxy is 
complex: co- and counter-rotating, as well as non-rotating, stellar 
orbits are needed simultaneously to reproduce the observed kin- 
ematics (Fig. |9). The NSC seems to be embedded in a disk that is 
most likely younger and more metal rich than the main body of the 
galaxy (Fig. lilt. Due to insufficient S/N we can only provide a con- 
servative upper limit for a possible black hole of 10 7 Mq. All these 
facts point to a complex formation history of the nuclear region 
in FCC 277. Most likely gas dissipation and merging played an im- 
portant role in shaping the nucleus of this galaxy. To check whether 
this is a common phenomenon among early-type galaxies, a larger 
sample is needed and can be obtained with current observing facil- 
ities. 
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